Novel Spin-Gap Behavior in Layered S = ^ Quantum Spin System TiOCl 
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We have investigated the spin fluctuations and their interplay with lattice instabilities in a layered 
S = \ quantum spin system TiOCl. Our '''^'''STi and "^^Cl NMR data demonstrate a pseudo spin- 
gap behavior below T* = 135if , followed by successive phase transitions at Td = 94A" and Tc2 = 
66 A" into the singlet ground state with an unusually large energy gap Eg — 430 A'. The broad 
distribution of the local electric field gradient in the intermediate phase between Tci and Tc2 suggests 
unconventional spin and lattice (probably orbital) states with an emerging spin-gap. 
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The provocative proposal that the mechanism of high 
temperature superconductivity in layered cuprates may 
be related to the exotic properties of low dimensional 
quantum spin systems, such as the RVB (Resonanting 
Valence Bond) state [1], has been a major driving force 
behind the rapid advance of quantum magnetism. Nat- 
urally, the initial emphasis was placed on understanding 
materials involving Cu^"*" ions with the 3d^ configuration 
{S = i) [2,3]. More recently, S = ^ quantum mag- 
nets involving [4-6] and Ti^+ [7-12] ions with Sd^ 
configurations have been attracting strong attention. Po- 
tential interests include the realization of doped metallic 
states in 3d^ Mott insulators, and ultimately, supercon- 
ductivity [6]. In addition, the near degeneracy of the 
t2g-orbitals often gives rise to the orbital and/or charge 
ordering [13]. In the orbital ordered state, the configu- 
ration of the occupied 3dxy,yz,z'x-0T:hitals exhibits a long 
range order. On the other hand, the valence of ions dif- 
fers site by site in charge ordered states. The additional 
orbital and charge degrees of freedom make the under- 
lying physics of 3d^ quantum spin systems more intrigu- 
ing, yet more complicated. A fascinating example is a 
mixed valence system NaV205 [5,14] (the average va- 
lence of V-ions is -1-4.5). NaV205 undergoes orbital and 
charge ordering at Tc = 34if , where a large energy gap 
Eg = 98K opens [14-17]. The ratio 2Eg/kBTc = 6 is 
much larger than the BCS value 3.5 expected for conven- 
tional spin-Peierls transitions with lattice-dimerization. 
Intensive theoretical and experimental efforts have been 
underway to understand the exotic phase transition. 

An equally exciting new avenue to investigate S — ^ 
quantum spin systems with 3d^ electrons is the titanates 
involving Ti^+ ions. Starting from the metal-insulator 
transition in Sri_a;La2,Ti03 [7], growing efforts are under 
way in search for a new form of quantum magnetism in 
titanates: strong orbital fluctuations in LaTiOa [8]; a 
quasi ID S" = i chain system (Na,Li)TiSi2 06 [9]; 3D 
Pyrochlore system MgTi204 [10]; and a quasi-2D layered 
system TiOCl [11,12]. 

Among these titanates, the TiOCl system has partic- 



ularly unique characteristics. First, Ti'^+O^" form bi- 
layers separated by Cl~ bi-layers. The quasi-2D layered 
structure may be considered a close analogue to those 
realized in high Tc cuprates. In fact, Beynon and Wilson 
[11] reported the very little temperature dependence in 
uniform susceptibility x{T), and discussed the possible 
realization of a RVB ground state in the Mott-insulator. 
They also noted that x(^) is very sensitive to impurities, 
and demonstrated that a Sc^+ ion {S = 0) substituted 
into a Ti'''+ site {S = \) gives rise to a localized spin 
S = i, in analogy with creation of a free spin by Zn^"*" 
{S — 0) substituted into Cu^+ [S — i) sites in cuprates. 
Second, in a very recent report, Seidel et al. observed 
a sharp, nearly isotropic drop of xC^) below ~ IGQK in 
defect-free samples, signaling the emergence of an energy 
gap [12]. The drop in xi^) is even more pronounced in 
high quahty single crystals. Since x(^) above ~ 
can be fitted with a 5 = ^ Heisenberg chain model with 
the nearest neighbor exchange interaction J — GGOif, Sei- 
del et al. proposed, based on LDA-I-U calculations, that 
the effective dimension of the TiO layers are reduced from 
2D to ID by an orbital order at Ti'^+ sites along the a- or 
b-axis. The gapped behavior of x{T) was attributed to 
a spin-Peierls transition [12]. Third, but related to the 
second point, the near degeneracy of the energy levels of 
different 3d^ orbital configurations without mixed-valence 
nature might make TiOCl an ideal model system to inves- 
tigate spin-Peierls- Zifce transitions with additional orbital 
degrees of freedom but probably without charge order. 

In this Letter, we report the first NMR investigation 
of TiOCl. While the spin-Peierls transition in CuGeOa 
{without orbital and charge degrees of freedom) and the 
exotic order in NaV205 {with orbital and charge degrees 
of freedom) have seen many detailed microscopic inves- 
tigations, to the best of our knowledge this is the first 
successful microscopic experiments reported for TiOCl. 
We demonstrate that TiOCl reveals a unique spin-gap 
behavior accompanied by lattice instabilities, and un- 
dergoes successive phase transitions at Tci = 94 ± 2K 
and Tc2 = 66 ± 2K. Unlike CuGeOg and NaV205, 
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the fluctuation effects in TiOCl are so strong above 
Tci that a pseudo spin-gap manifests itself as high as 
T* = 135 ± lOK. Moreover, our observation of a broad 
continuum in the NMR Hneshape data indicates that the 
intermediate phase between Td and Tc2 is not a simple, 
dimerized state in ID. Below Tc2, TiOCl undergoes a 
first order phase transition to a fully gapped state (rather 
than to the gapless RVB state) with an extraordinarily 
large energy gap Eg = 430 ± 60K. The unusually large 
energy gap {2Eg/kBTci.c2 = 10 ~ 15) as well as the 
presence of the intermediate state between T^i and Tc2 
suggest that the observed phase transitions are not con- 
ventional spin-Peierls transitions, and point towards the 
significance of the roles played by the additional orbital 
degrees of freedom over the entire temperature range. 

Our TiOCl single crystals were synthesized by stan- 
dard vapor-transport techniques from Ti02 and TiCla 
[12]. A large number of very thin, flaky single crystals 
with typical dimensions of 2mm x 2mm were assembled 
on a Macor sample holder. NMR measurements were 
conducted by applying an external magnetic fleld in par- 
allel with the aligned crystal c-axis. Random alignment 
within the TiO-plane prevented us from conducting mea- 
surements along the a- and b-directions. We empha- 
size that, unlike ^^'^^Cu and ^^V NMR, the sensitivity 
of 47,49fpj NMR is notoriously low [18]. Very short trans- 
verse relaxation times T2 at Ti sites and the small volume 
of the available crystals made the NMR measurements 
even more difficult. We needed to average spin-echo sig- 
nals up to 10^ scans to obtain a reasonable signal to 
noise ratio. In general, NMR signal intensities increase as 
1/T in proportion to the Boltzman factor. However, the 
exponentially growing spin-lattice relaxation time Ti be- 
low Tci slows down the NMR pulse-sequence, hence the 
measurements at lower temperatures were equally diffi- 
cult and time consuming. Despite our intensive efforts, 
we have been able to find 47,49rjij ^MR signals only for 
the central transition from the = +^ to — ^ state. 
The detection of ^^Cl NMR signal was somewhat easier, 
and we did manage to find the 1^ = ±| to ±i transi- 
tions above T^i [19]. All the nuclear spin-lattice relax- 
ation data were deduced by fitting the nuclear spin recov- 
ery after an inversion 7r-pulse to the standard rate equa- 
tions. We confirmed that ^^Ti (nuclear spin / — 5/2) and 
■^^Ti (nuclear spin / = 7/2) NMR gives identical 1/Ti. 
Given that the magnetic recovery process of nuclear mag- 
netization M{t) is dominated by different terms for ^^Ti 
{M{t) ~ exp[-15t/Ti]) and^OTi {M(t) ~ exp[-28t/Ti]), 
we conclude that 1/Ti is dominated entirely by magnetic 
fluctuations at all temperatures. The difference in the re- 
covery characteristic of M{t) also helped us identify the 
pairs of ^^Ti and ^^Ti lines below T^a, where the NMR 
lines split into doublets. 

The experimental information on the spin degrees of 
freedom is summarized in Fig. 1. Quite generally, the 
47,49fj'j nuclear spin-lattice relaxation rate, caused by low 



frequency magnetic fluctuations may be expressed as [20], 

where F{q) is the form factor of the electron-nucleus hy- 

perfine couplings, and x"(q, ti;„) is the imaginary part of 
the dynamical susceptibility at the observed NMR fre- 
quency iVn- ^/T\ appears to asymptote to a constant 
value near 300i^, l/Ti 1400 sec'^ This imphes that 
x"(q, Wfi) asymptotes to a Curie law at higher tempera- 
tures. This is consistent with the behavior expected for 
typical S' = i ID Heisenberg chains such as SraCuOa 
[21]. However, we caution that l/Ti ~ const, may also 
be expected iov S = \ 2D Heisenberg model [22]. Cer- 
tainly x(T) of TiOCl also fits nicely to the ID Heisenberg 
chain model between and SOOif [12], but we recall 

that NaVaOs turns out to be a 2D spin-charge-orbital 
hybrid system [23] despite the similarly nice fit of x(^) 
to the ID model. 

In general, the growth of short-range order enhances 
low frequency spin fiuctuations (hence l/TiT) with de- 
creasing temperature in the absence of a gap in the spin 
excitation spectrum. The most striking feature in Fig.l 
is that 1/TiT begins to decrease below T* ~ 135-?^. 
This implies that low frequency spin fluctuations are sup- 
pressed below T* , by almost two orders of magnitude 
between 135-?^ and Q'bK. Interestingly, the observed be- 
havior of 1/TiT is qualitatively similar to the pseudo-gap 
phase in underdoped high cuprates. Needless to say, 
it does not necessarily imply that the underlying mecha- 
nism is identical. 

We can also probe the lattice degrees of freedom by 
observing the EFG (Electric Field Gradient) reflected on 
NMR lineshapes. In Fig. 2, we present the 47,49rpj ^-^^^ 
^^G\ NMR lineshapes for the central transition. We ob- 
served a single NMR line for both Ti and CI down to 
Tci = 94 ± 2K. For '^"''Cl, we also managed to detect 
Iz = ±1 to = ±j satellite transitions, but again 
there was only one kind of signal [19]. These results in- 
dicate that there is only one kind of Ti and CI site in 
TiOCl within our experimental resolution, and rule out 
any potential orbital order configurations above; T,.i that 
would lead to more than one inequivalent sites. Below 
Tci, both 47,49Ti and ^^Cl NMR central fines begin to 
broaden, signaling a second order phase transition. We 
confirmed that the mag nitude of the ^^Cl NMR fine split- 
ting is inversely proportional to the external magnetic 
field. Therefore the line splitting is caused by the second- 
order nuclear quadrupole interaction with the EFG. We 
also found that the drop of 1/TiT is accelerated below 
Tel. In addition, close inspection of the x(T) data re- 
ported in [12] reveals a kink at T^i , followed by a rapid de- 
crease with temperature with positive curvatures. These 
results suggest that a clear gap structure emerges in the 
spin excitation spectrum at Td, accompanied by a static 
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distortion of the lattice. The NMR hneshape exhibits 
a broad continuum below Td down to Tc2 as summa- 
rized in Fig.3(b). The broad continuous distribution of 
the local EFG environment implies the presence of nu- 
merous inequivalent Ti and CI sites in the TiOCl lattice. 
The most plausible scenario is that the emergence of the 
spin-gap at Td is accompanied by an orbital-order (pos- 
sibly incommensurate). That is, the phase transition at 
Tci is not a simple spin-Peierls transition due to lattice 
dimerization, even though the well-developed short-rage 
spin order below 400K and the gapped behavior of x{T) 
below ~ lOOK [12] suggest the contrary. 

We found additional evidence for the involvement of 
the lattice in the phase transition at Td by ^^Cl nu- 
clear relaxation data as shown in Fig.3(a). ^^1/Ti 
is two orders of mag nitude slower than 1/Ti at 47,49rpj_ 
sites. This implies that the magnetic hyperfine coupling 
of ^^Cl nuclear spins with Ti 3d^ electron spins is at 
least by an order of magnitude smaller. The qualita- 
tively different temperature dependence between "^^l/Ti 
[see Fig.3(a)] and 1/Ti [see Fig.l(a)] suggests that slow- 
ing of EFG fluctuations below T* toward the phase tran- 
sition at Td causes the cusp of ^^1/Ti at Td [24]. Even 
more straightforward evidence is that the ^^Cl spin-echo 
decay (transverse T2 relaxation) in Fig. 3(c) and (d) shows 
typical motional narrowing effects below T* towards Td 
due to lattice softening. The slow lattice dynamics at 
or below the ^^Cl NMR frequency (37.6 MHz) trans- 
forms the spin echo-decay to Lorentzian-like as Td is 
approached. Taken together with the precursive sup- 
pression observed for 1/TiT below T* in Fig. 1(a), we 
conclude that the drastic softening of the lattice below 
T* drives the observed pseudo spin-gap behavior prior to 
the actual second-order phase transition at Td = 94if . 
The high onset temperature ^ 200K of the gradual de- 
crease of 1/Ti suggests that phonon softening at higher 
frequency scales begins at much higher temperature. 

At Tc2, the continuum of the ^^Cl and 47.49 j^j^lR 
lineshape suddenly collapses into doublets, whose line po- 
sitions correspond to the two extrema observed above 
Tc2- This strongly suggests that the imit cell of Ti- 
OCl has two and only two inequivalent Ti and CI sites 
below Tc2 [25]. The low- frequency Ti 3d^ spin fluc- 
tuations exhibit an activation behavior below Tc2, as 
shown in Fig. 1(b). From the fit to an exponential form, 
l/TiT ^ exp(—Eg/kBT), we deduce the energy gap 
Eg = 430 ± 60K. These results are consistent with 
the formation of a singlet ground state at Tc2 by lat- 
tice dimerization. Our results are also consistent with 
the recent observation of doubling of the unit cell along 
the b-axis at 66 ± IK by Y.S. Lee et al. based on x- 
ray scattering measurements [26]. On the other hand, 
there are many experimental signatures which are at 
odds with conventional second order spin-Peierls tran- 
sitions. 1/Ti [Fig.l], the NMR lineshapes [Fig.2], and 
35l/ri [Fig.3(a)] change discontinuously at Tc2. ^^Cl 



spin-echo decay curve also changes suddenly to a slow 
Gaussian type at Tc2 [Fig. 3(d)]. This indicates that the 
soft vibration of the lattice disappears suddenly below 
Tc2. AH these NMR results as well as the discontinu- 
ous change in x{T) [12] and the x-ray data [26] at Tc2 
indicate that the phase transition at Tc2 is first order. 
Furthermore, the observed energy gap is unusually large, 
2Eg/kBTd,c2 = 10 ~ 15. In the case of conventional, 
mean-field spin-Peierls gap, one would expect the ratio to 
be close to the BCS value, iEg/ksTc = 3.5, as observed 
for CuGeOg {2Eg/kBTc = 3.3 [27,28]). The observed 
ratio in TiOCl is even larger than that for NaV205, 
2Eg/kBTc ^ 6 [14,17]. It is important to reahze that Eg 
is comparable to the apparent exchange energy J = 660iir 
estimated by the ID fit of x(^) [12]. Such a large mag- 
nitude of the energy gap Eg strongly suggests that the 
spin excitations from the singlet ground state is dressed 
by other electronic degrees of freedom, most likely of the 
orbital origin. 

To summarize, TiOCl is a rather unique layered S = ^ 
quantum-material with a pre-existing pseudo spin-gap 
above Td , the unconventional intermediate spin and lat- 
tice (probably orbital) states between Td and Tc2, and 
a first order phase transition into a singlet ground state 
with an unusually large energy gap. At first sight, the 
strong fluctuation effects evidenced in the pseudo spin- 
gap behavior suggest a genuinely ID nature of the TiOCl 
along the a- or b-axis, achieved by an orbital order well 
above Td- We recall that the pre-existing orbital order 
above Td may reduce the effective dimensions of the TiO 
layer from 2D to ID, as suggested by Seidel et al. based 
on LDA+U calculations [12]. On the other hand, the 
additional orbital degrees of freedom along the two or- 
thogonal directions within the 2D TiO-layers may cause 
strong orbital fluctuations, which could effectively sup- 
press the tendency towards a spin-Peierls transition at 
finite temperature. In such a scenario, the manifestation 
of the pseudo spin-gap above Td may be the consequence 
of the suppression of the spin-Peierls transition, and the 
successive phase transitions at Td and Tc2 may arise from 
a competition between different orbital orders. In addi- 
tion, if the second nearest-neighbor exchange interaction 
in the zig-zag chain structure along the a-axis (see 
Fig. 1(c) in [12]) exceeds 0.241Ja (where Ja being the 
nearest-neighbor exchange interaction along the a-axis), 
a spontaneous spin dimerization along the a-axis may be 
favoured [29] over the unit-cell doubling along the b-axis 
observed below Tc2, introducing frustration. In any case, 
the unprecedented behavior of TiOCl points towards the 
crucial roles played by the orbital degrees of freedom and 
their fluctuations below ~ 200K. Recalling that TiOCl 
is not a mixed- valence system such as NaV205, TiOCl 
may be an ideal model spin-Peierls system with addi- 
tional orbital degrees of freedom, and we call for futher 
microscopic studies. 
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FIG. 2. NMR lineshapes of the = +i to -i central 
transition for ^^Cl (left panels) and (j-^gj^t panels). The 

FFT envelope was obtained at 9 Tesla by repeating spin-echo 
measurements at several frequencies. 
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FIG. 3. (a) ^^l/Ti. (b) ^^Cl NMR frequencies observed 
at 9 Tesla. Dashed lines represent continuum. Solid er- 
ror bars represent the frequency range corresponding to the 
half-intensity, (c) *^C1 NMR spin-echo decay observed at the 
main peak, as a function of the pulse separation time 2r be- 
tween 90-180 degree R.F. pulses above Td, and (d) below Td. 
Solid curves are guides for eyes. 



FIG. 1. (a) 1/Ti (o) and 1/TiT (•) at '''^'^Ti sites. 1/Ti 

was measured at the center peak for T < Td- Solid and 
dashed curves are guides for eyes, (b) The same 1/T\T data 
plotted in a semi-log scale. Solid line is the best exponential 
fit. 1/T\T was identical for the doublets below Tc2. 
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